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UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER LABORATORY 
STATISTICAL LIBRARY 

KSL 2.03 - 513 
Means J standard, deviations^ third and fourth moments about 
the means 
Entire program. 

V < 1^5 -where v is the number of variables)!: no practical 
limit on q^ the sample size 

For each of a set of v variables^ this routine will calculate 
the mean, the standard deviation, and the third and fourth 
moments about the mean. Much can be learned about the sample 
distribution from these four statistics. The mean is a mea- 
sure of central tendency. The standard deviation is a measure 
of the dispersion of the distribution. The third moment 
indicates the amount and direction of skewness » The fourth 
moment is a measure of the kurtosis of the distribution = 
The routine will read data either in the form of signed 
fractions or in the form of unsigned single digits (O, 1, ... 9)» 
If means and standard deviations only are desired, computer 
time can be saved by suppressing the calculation of the third 
and fourth moments. In any event, this routine is preferred 
over K-I7 (also K-8, KSL 2.0l}for'the v(v + l)/2 cross-products 
are not calculated o For a large v, the saving in computer 
time may exceed. 90^» 

lo Read, master program. Routine stops on ^kodk^ A stop on 
FFOO3 indicates a sum check' failure. 

2, Read parameter tape. Routine stops on S4o4n. 

3. Read data tapeo Routine ends on 2^084. 

Another problem can be begun by reading a new parameter tape. 
PARAMETER TAPE PREPARATION: 

There are four parameters (r, d, q, and p) for each problem. 
The parameters are punched on tape with sexadecimal terminating 
symbols as follows: 
r N d-" J q F p L. 



METHOD OF USE: 
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To read data consisting of single unsigned digits, set r = 0. 

fo read data consisting of signed fractions, set r = 1. 

Set d equal to the number of decimal places desired in the 

means and standard deviations. The third and fourth moments 

are always punched to 10 decimal places . 

q is the sample size or number of rows of data. 

To suppress the third and fourth moments, set p = 0. To 

calculate and punch the moments, set p = lo 

DATA TAPE PREPAEATION; 

The data tape is punched by rows and consists of q rows 

of V variables each. Each row of the data tape is terminated 

by an N symbol. If an F follows a row instead of an N, the 

computer will stop and another section of the data tape can 

be inserted in the reader. By raising the black switch, the 

reading of the data tape is resumed. 

When r =0, an element of a row must be a single unsigned 

digit (O, 1, ... 9)« If these are considered as integers, 

then the scaling on each will be 10" . 

When r = 1, each row element must be punched as a signed 

fraction with any number from zero through 12 decimal digits. 

SCALING IN THE RESULTS: 

If the scaling on variable j is 10"-^j, then in the results 
the means and standard deviations will also bp- scaled by 
10 -^j . The third moment will be scaled by 10~ -^j and the 
fourth moment will be scaled by 10 ^j . 

THE FORM OF THE RESULTS? 

The means and standard deviations are printed out in parallel 
columns tenninated by an N. 

If moments are also calculated, these will follow in a second 
set of parallel columns terminated by an N. 

A THREE-VARIABLE EXAMPLE: 

Parameter tape 
ON3J6f1L 
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DURATION; 



FOEMUIAS: 



Data Tape 

232N 5^2N 15IN 06ON ^TON 371N 

Results 

^+217 +13^ 
+533 +1^9 
+100 +082 
N 



X.'s 



s . s 

J 



M.5's 



r^-ooojkojko^ +0006025^62 *! 
I -0009259259 +000807^07^ f 
I +0000000000 +0000666666 / 



M.^'s 



^K 



sec,/d = decimals | 

q = samples I 

^= variables y 



Read master tape 
Read data in 
Calculation 
Pimch 



th 
Let X. . be the i observation (i =1, 2. 



30 

.ooi^(d + i)-|rq. 

.003'^q 

.0357^((i + 1) means and standard 

deviations 

,k^y^rdL and i|-th mom 

N , , . th 
. q; on the j 

variable (j =1, 2, ... v). All summations below are over 

i from 1 through q. 

Mean (central tendency) 

X. = ^ Z X. . 

standard deviation (dispersion) 

s. = (M.^)^/^ =l^z\x.. - X.)^ 

Third moment {skevness) 

M.^ =iE (X.. -x.)5 

Fourth moment (kurtosis) 

M,^=iE(X.. -x.)*^ 



1/2 



- k , 



3 
If the distribution is symmetric, M. = 0» 

When M. is negative, the distribution is skewed to the left: 

when positive, the distribtuion is skewed to the righto For 

purposes of comparison, iise 

's 



A. = M.^r ^ 



For a flat distribution, M. will tend to be large; for a 

k ^ 
steep distribution, M. will tend to be small. For purposes 

of comparison, use B, = M. /s . « 

J J J 

For a normal curve, A. will be equal to zero and B. will be 

J J 

equal to three , 

1« The routine compares the number of variables in subsequent 

rows with the number in the first row. If these do not agree 

for any row, the computer will stop on FFOOO at location OSl. 

2 
2. If the variance, s. , is negative due to rounding errors 

and outside the tolerance limit (lO" ), the machine will stop. 

on FF0S6 at location 0N6» If a negative variance is vithin 

the tolerance limit, it is set to zero and the calculation is 

continuedo 



DAT E October 12, I96Q 

PROGRAMMED B Y Syx.£.^cLcL^ y-Xo^zM ^ 

APPROVED B Y lyyl JUiAA/d<A J 




nj 



LOCATIOI^ 


ORDER 


NOTES PAGE 1 KSL 2.03 j 


Abs. 


Rel.lsym. ! 


00 UK 






1 

1 


11 


(2) 


00F;002F 








(10) 


OOF OOIOF 








(1^5) 


00145F 00145F 








(299) 


FF299F 002 99F 








(CL) 


50F itll02ij-F 




16 




(J) 

(D2) 

(d6) 
(z) 


OOF OOOOOOOOOOOIOOJ 
OOIF OOF 
6F299F 002 99F 
OOF OOF 




20 


9 (Pl6) 


OOK 




7^ 

1^ 




(N12 ) 


OOK 




115 




(Rl) 


OOK 
OOK 


Read in digits 


12i^ 





(N13) 


K5F i|-27L 
h6^L 50 (Z) 
QlkF LO(lO) 
32 7L Li|(lO) 
66(10) S5F 




129 






i|-OF?L55L 
iM- (D2 ) 465L 




151 


7 




221L 22F 
OOK 




132 





(L) 
(LI) 


193F 401F 

92575F ii-lF 
8i4f lo(io) 

32(12) lA(lO) 
50F 7^(10) 
S5F 4of 


Read parameters 


138 




(L2) 


26(Ll) i^2l(L2) 
L5F 40F 
L51F L^IF 




141 

i- 


9 

1 ,... , . 




401F 321 (L) 





LOCATION 


ORDER 


NOTES PAGE 2 KSL 2 .03 


Abs,|Rel. 


Sym. 


l44f 4o6f 




l42 jio 










50 (z) 1938F 










664f S5F 










4o8f 50 (z) 
L55F 0039F 
405F L52F 
0039F 402F 
L53F 0020F 
46(?i) 46(P2) 
92139F L52f:^ 




152 


20 




36(L3) 26(l4) 








(L3) 


50299F 50 (L3) 
24(N13) L0(2) 
L0(2) 343 (L3) 
L55(N13) 263(14) 


Read first raw digits 

■ 






(14) 


50299F 50 (l4) 
24(N12) L0(2) 
302 (l4) L52l(N12) 
409F 1020F 
L0(299) 4210F 


Read fractioAs 
NumlDer of varia'bles 


162 


30 


(L5) 
(l6) 


L5(299) l4(i45) 
42 (L5) 92707F 
JOF 4lF 
■f5(L5) 42 (L5) 
LO(CL) 32 (L5) 
4l7F 26 (l8) 
50299F 50 (l6) 
26(10.3) L0(2) 
K)(2) 343(l6) 
L55(N13) 263 (L7) 




172 


40 


(L7) 


50299F 50 (L7) 
26(N12) L0(2) 
302 (L7) L521(KL2) 




175 

V 


43 




L09F 40F 





LOCATION. 


ORDER 


NOTES PAGE 3 


KSL 2.05 


AbSo 


Eel. 


Sym. 


L3F 36 (l8) 






176 


44 








(T) 


FFF OOF 


incorrect number VHriables R 


.H.A. has tally 






(l8) 


L5(299) 40 (T) 
42(]^6) L5(T) 
46 (Lll) 463(111) 












464(1.11) 469(111) 






182 


50 


(Lll) 


l4(i45) 421(111) 
462 (Lll) l4(i45) 
424 (Lll) 465 (Lll) 
l4(i45) 427 (Lll) 
468 (Lll) l4(i45) 
429 (Lll) 46 (L9) 
50F 7J8f 
40F l4f 
40F 50F 
7JF 40F 






192 


60 




50F l4f 
40F L55F 
32 (L9) IJF 

40F l4f 
4qf 50F 
7JF l4f 


2 
Store E 2^ = 1^2 

Store Z - = r^ 
4 








(L9) 


40F F5(t) 
l4(D2) 4o(t) 
46(d6) lo(d6) 

LOIOF 321 (l8) 


Store E 1 =y-^ 




202 


70 


(110) 


F57F 42 7F 

lo6f 322 (lio) 

26 (S) L52F 
36 (l6) 26 (L7) 

OOK 






206 





(-) 


L5(299) L4(i45) 
4o(t) 42 (d6) 


3 4 
Compute M. and M, 

J 3 




208 

i 

i — »■ 


2 




L5(t) 46i(s4) 







L0C4TI0N, 


ORDER 


1 

NOTES PAGE 4 


KSL 2.03 


AbSa 

209 


Rel. 


Synie 


465 (S4) i^25(s4) 






5 










i^215(Sil-) i^2l6(s4) 












IA(li^-5) i^-2i^(s4) 












i^6lif(si^) 46 (Sl) 












42 (S2) 421 (S5) 












L4(i45) 421 (S4) 












4218 (S4) 4619 (s4) 






2X6 


10 




l4(i45) 4211 (s4) 
4615 (S4) 26 (s4) 










(sif) 


L55F 56(S1) 
50F 75F 
0020F 402F 
5OF 7JF 
40F 5OF 
75F OOIF 
401F OOIF 
l4if 405F 


M^^ =7^- 4x1^-+ 6x^X2 ' 


3^ 


226 


20 




50F 75F 
401F OOIF 
l4if 4o4f 
L55F l4f 
L02F L04f 

4oF 50F 
75F 4oiF 
OOIF 

l4if 

402F 5OF 
75F OOIF 


Store M^ 

M^^ =-^ - 5x1^ + 2x^ 




236 


30 




L02F L4f 
40F 26 (Si) 


Store M.^ 








(SI) 


L5F LOF 
401F 4lF 

L5rF 52 (S5) 


(7-2-5?) 




2ia^ 


55 




L4(j) 52(S2) 









i^OCATION 



OEDER 



NOTES 



PAGE ^ 



KSL 2.05 



Abs 



Rel. 



Sym. 



2i^2 



36 



246 
248 

2^9 



40 
42 





259 



269 



274 



10 



20 



25 



(S2) 
(S3) 



(P) 



(PI) 



(P2) 



(PR) 



(P3) 



(pM 



FF2F 4lF 
262(33) 50(33) 

26 (ri) 4pF 

F5(t) L4(D2) 

4o(t) 46 (d6) 

L0(d6) LOIOF 
362(3) L5(299) 

OOK 

l4(i45) 4o(t) 
42 (d6) L5(t) 
423 (P) l4(i45) 

422 (pi) L5F 
5012F 50 (Pl) 
26(pi6) 92963F 
JOF L5F 
5012F 50 (P2) 

26(Pl6) 9213IF 
92515F F5(t) 
l4(d2) 4o(t) 
46 (d6) lo(d6) 

LOIOF 321 (P) 
9277OF 92535F 
L55F 32 (P5) 
92147F L5(299) 
l4(i45) l4(i45) 
l4(i45) 4o(t) 
42 (d6) L5(t) 
425 (PR) l4(i45) 
422 (P3) L5F 
501OF 50 (P3) 
26(pi6) 92965F 

JOF L5F 

5010F 50 (p4) 
26(Pl6) 9213IF 



store s 



Print routine 



LOCAiriON 



ORDER 



NOTES 



PAGE 6 



KSL 2.03 



R@l. 




279 


30 


280 


31 


281 





285 


2- 



Sym. 



(P5) 



9^515F F5(t) 
IA(D2) il-O(T) 
46(d6) L0(d6) 
LOIOF 323 (PR) 
9277OF 92575F 
9213IF feil-(L) 

OOK 

L3F 5^(L) 
FF3F 26 (l) 
F72889F 211040F 
26L 26IN 



Sum check 



